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A proper treatment of the resonant continuum is to take account of not only 
the energy of the resonant state, but also its width. The effect of resonant states 
on pairing correlations is presented in the framework of the relativistic mean field 
theory plus Bardeen-Cooper-Schrieffer (BCS) approximation with a constant pairing 

o 

■ strength. The study is performed in an effective Lagrangian with the parameter set 

o' 

C*~> ■ NL3 for neutron rich even-even Ni isotopes. Results show that the contribution of 

O 

the proper treatment of the resonant continuum to pairing correlations for those 
nuclei close to neutron drip line is important. The pairing gaps, Fermi energies, 

<■> : 

pairing correlation energies, and binding energies are considerably affected with a 
proper consideration of the width of resonant states. The problem of unphysical 

•i-H . 

particle gas, which may appear in the calculation of the traditional mean field plus 
BCS method for nuclei in the vicinity of drip line could be well overcome when the 
pairing correlation is performed by using the resonant states instead of the discretized 
states in the continuum. 
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I. INTRODUCTION 

The study of properties of exotic nuclei, in particular the structure of nuclei far from the 
/3-stability line, has recently attracted wide interest both experimentally and theoretically. 
In those nuclei, the presence of nucleons, which have separation energies appreciably smaller 
than those in /3-stable nuclei and have a far extending nucleon distribution, leads to various 
interesting physical phenomena. Due to the closeness of the Fermi surface to the particle 
continuum in exotic nuclei the coupling between bound states and the continuum become 
important. The description of exotic nuclei has to explicitly include the coupling between 
bound states and the continuum. There exist several extensive studies of exotic nuclei in 
relativistic and non-relativistic microscopic approaches 0, 0, 0, |^. A key ingredient of those 
models is how to properly treat the pairing correlations which have an important influence 
on physical properties in exotic nuclei. In general, pairing correlations in open shell nuclei 
can be treated by the Bardeen-Cooper-Schrieffer (BCS) theory or through the Bogoliubov 
transformation The main feature of the Bogoliubov transformation, compared with the 
simple BCS theory, is that the Hartree-Fock equation or Dirac equation and the gap equa- 
tions are solved simultaneously with self-consistent fields. Both of them could give a good 
description of the pairing correlation if the nucleus is not too far away from the /3-stable 
line|6l. The simple BCS method mav not be reliable near the drip line because the continu- 

riQ 

ous states were not correctly treated |6|, LD] - The contribution of the coupling to the continuum 
would be prominent when the nucleus close to the drip line, therefore a proper treatment of 
the continuum becomes more important. Single particle resonant states in the continuum, 
which differ from bound states, are metastable states with life times due to a sufficiently 
high centrifugal barrier (for neutrons and protons) and Coulomb barrier (for protons). The 
effect of the widths of resonant states on pairing correlations was not estimated by the usual 
non-relativistic Hartree-Fock-Bogoliubov (HFB)^,!^ or the relativistic Hartree-Bogoliubov 
(RHB) 81. The effect of a proper treatment of the resonant continuum on the pairing in the 

nnn 

non-relativistic approach has been investigated [9 1 . llflL II 1| . The investigation in the relativis- 



tic approach is reported very recently 
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In this paper, we aim at the investigation on the effect of the resonant continuum on 
pairing correlations of neutron rich nuclei in the relativistic mean field theory plus BCS 
(RMF+BCS) approximation. Particularly we shall focus our attention on the effect of 
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widths of single particle resonant states in the description of neutron rich nuclei. It has 
been pointed out that pairing correlations could be well described by the simple BCS theory 
if single particle states in the continuum are properly treated^, Q, 11 ] . It is known that 
narrow resonant states in the continuum are well localized inside the nucleus due to a 
sufficiently high centrifugal barrier and Coulomb barrier (for proton only). In our previous 
work[13], we have shown that the nuclear dynamical processes can also be well described by 
the simple method where the continuum was replaced by only a few single particle resonant 
states. Therefore, we may expect that using a few single particle resonant states instead of 
the discretized continuum, the pairing space would be large enough in order to give a good 
description on the ground state properties of neutron rich nuclei in the RMF, even within 
the BCS approximation. 

Usually pairing correlations in the BCS approximation are treated in a phenomenological 
way based on empirical pairing gaps deduced from odd-even mass differences. This procedure 
works well in the valley of /3-stability, where experimental masses are known. The predictive 
power for pairing gaps for nuclei far from the /3-stability line is thus limited due to the un- 
known masses of such nuclei. Therefore, a constant pairing strength is usually employed in 
the BCS calculation [14]. We shall investigate the pairing effect in the RMF with the param- 
eter set NLs|l^) which give a good description of not only the ground state properties jl^] 



but also the collective giant resonance 
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201 ] . To investigate the width effect of sin- 



gle particle resonant states on pairing correlations, we compare the results in three types of 
calculations: (a) the RMF+BCS approach where quasi-particle wave functions are obtained 
by a box boundary condition; (b) the extended RMF+BCS approximation where wave func- 
tions of resonant states are calculated by imposing a proper scattering boundary condition 
for the continuous spectrum but the width of resonant states is not considered, which is de- 
noted as RMF+BCSR; (c) the full RMF+BCS (RMF+BCSRW) calculation where widths 
of single particle resonant states are taken into account explicitly. 

The paper is arranged as the follows. In Sec. II and Sec. Ill the relativistic mean field 
theory and the BCS approximation with explicitly considering the widths of the resonant 
continuum are presented. The effect of the resonant continuum on pairing correlations is 
formulated by introducing a continuum level density into the gap equations. The results 
calculated by different models mentioned above for neutron rich even-even Ni isotopes are 
given and discussed in Sec. IV. Finally we give a brief summary in Sec.V. 
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II. RELATIVISTIC MEAN FIELD THEORY 



The RMF theory 2l|, |22( is based on the effective Lagrangian density which includes 
the nucleon field (ip), the isoscalar scalar meson field (a), the isoscalar vector meson field 
(uj), the isovector vector meson field (p), and electromagnetic field (A). The Lagrangian 
density is expressed as the following: 

-g^ail) - g^ip^u^ip - gpip^pf, ■ t4> 

-e^A, i> , (1) 
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where M, m CT , m u , and m p denote the nucleon, cx, cu, and p mesons masses, respectively, 
while g a , g u and g p are the corresponding coupling constants for the mesons, respectively. 
The vectors in isospin space are denoted by bold-faced symbols. , R tiv and F^ v are the 
field tensors of the vector fields ui, p and of the photon, respectively, they are defined as: 

QT = d»uj v - d u u» , (2) 

R l» = Q» p v _ QV p » ? (3) 

In the Eq.(l), a non-linear scalar self-interaction termj^S] U(cr) of the a meson has been 
taken into account: 

U{a) = \mlo* + l -g 2 o* + \g^ , (5) 

The Dirac equation of a single-particle state in the RMF approximation can be expressed 
as following: 

[ a . p + V (r) +/3(M- S(r))] tfj a = E a ip a , (6) 

where E a and ip a are the single-particle energy and wave function, respectively. S(r) and 
V(r) are an attractive scalar and repulsive vector potential, respectively. In the RMF ap- 
proximation, scalar and vector potentials are produced by the classical meson fields: isoscalar 
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a, 00 and isovector p mesons as well as the photon, which are obtained in a self-consistent 
calculation for the nuclear pound statefl fl Q. The nucleon spinor ip a in the Dirac 



equation is expressed as: 



1 / iGJr) 

Mr) = - S«,m(f)Xi , (7) 



r 



Fair) a ■ f 



where G a and F a are upper and lower components of nucleon wave function, \\ is the 
isospinor and $ K , m is a spherical harmonics function, a is a set of quantum numbers a = 
in, l,j, m, r 3 ) = (a, m, r 3 ). k is the Dirac quantum number given by: 

-(j + i), for J = l + \ _ 
+ (j + i), for j = I - \ 

For spherical nuclei, the Dirac equation can be reduced to coupled equations of the radial 
part for G a {r) and F a (r): 



d 


' G a (r) \ 




dr 


\ Fair) j 


H 



(9) 



-* M + E a -S(r)-V(r)) \ I G a {r) 

M - E a - Sir) + V(r) & ) \F a (r) 

where S(r) is the scalar potential: 

S(r) = g a a(r) , (10) 

and V(r) denotes the vector potential: 

V(r) = gjJoif) + g P r 3 poir) + e 1 n A Q (r) , (11) 

The meson and electromagnetic fields obey the radial Laplace equations: 
d 2 2d 

-7^cr(r) + -— <r(r) - m^cr(r) = -# CT p s (r) - g 2 (J 2 (r) - g 3 a z {r) , (12) 
dr A r dr 

d 2 2d 

^wo(r) + -—uj ir) - m^ (r) = -g^pvir) , (13) 
cp 2d 

^Po(r) + —^Po(r) - m 2 p p (r) = -g p p 3 (r) , (14) 
J^o(r) + l^Mr) = -ep c (r) , (15) 
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with 



= E 2j ^(\Ga(r)f - |F a (r)| 2 ) , (16) 

a 

*(r) = E -^(I^WI 2 + l^(r)| 2 ) , (17) 

a 

occ „ . 1 

= E -^^(|G a (r)| 2 + |F a (r)| 2 ) , (18) 

a 

= E (|Ga(r)|2 + |Fa(r)|2) ' (19) 

The sums a run over all occupied states. The Dirac equation Eq. (9) and the correspond- 
ing meson field functions Eq.(12 - 14) as well as the photon field function Eq.(15) with 



expressions for source terms (16 - 19) can be solved sel: 



relativistic mean-field theory can be found in Refs. 



ii, 



21 



consistently. More details on the 
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III. BCS APPROXIMATION WITH RESONANT CONTINUUM 

It is well known that the pairing correlation plays an important role in describing t 
ground state properties of open shell nuclei. The non- relativistic HFB0, 0] or RHB 8j] 
theory can provide a unified description on the mean field and pairing correlations, but they 
require much more numerical effort than the simple BCS theory. However, the width effect 
of resonant states in the gap equation is missing in most previous works of the BCS theory 
or the Bogoliubov transformation. Recently, properly treating the contribution of resonant 
continuum to pairingcorrelations has been investigated in the non-relativistic HF+BCS as 



.^correlations nas been investigated m tne 
-13, IlD, U, 24]. It has been pointed outQ 



well as HFB in Refs.|9j, lld,lll|, |24(. It has been pointed out|24( that the resonant continuum 
HF+BCS approximation with a proper treatment of resonant continuum can reproduce the 
pairing correlation energies predicted by the continuum HFB approach up to the drip line. 
In this paper, we investigate the width effect of single particle resonant states on the pairing 
correlations in the relativistic approach. The pairing correlations of open shell nuclei are 
treated based on the RMF within the BCS approximation. A constant pairing strength 
is adopted in the BCS theory. Although the extended RMF+BCS theory formulated with 
correct boundary conditions for the continuous spectrum is available in the literature [l^. 
in the following we briefly summarize the essential points together with our model. From 



Ref. jfj, assuming constant pairing matrix elements in the vicinity of Fermi level, one can 
obtain the so-called gap equations: 



1 



v/(e a -A) 2 + A 2 G 



iV 



(20) 



(21) 



v/(5 a -A)2 + A 2 _ 

where A, A, and G represent the Fermi energy, pairing gap, and the pairing force constant, 
respectively. N is the number of neutrons or protons involved in the pairing correlations. 
The solution of those two coupled equations, Eq.(20) and (21) allows one to find A and A. 
The contribution of the pairing interaction to the total energy can be expressed as following: 



E 



pair 



-G 



-i 2 



(22) 



where v 2 is the occupation probability of a state with quantum numbers a = (n, I, j, m, r 3 ) = 
(a, m, T3) and u 2 a = 1 — v \. For spherical nuclei, the eigenvalues of Dirac equations are of 
degeneracy for magnetic quantum numbers m within the same (n,l,j). One has: 



£a - A 



V(e a -A) 2 + A2 



(23) 



In principle, the continuous spectrum is also the solution of Dirac Eq.(9), both continuous 
and bound states constitute a complete set of basis. In most previous theoretical nuclear 
structure calculations, single particle states in the continuum are usually treated in a dis- 
cretization procedure by expanding wave functions in a set of harmonic oscillator basis or 
setting a box. This approximation, however, can be justified for very narrow resonances and 
gives a global description of the contributions from the continuum. In this work we intro- 
duce single particle resonant states into the pairing gap equations instead of the discretized 
continuous states in order to investigate the width effect of resonances on pairing correla- 
tions. The resonant state wave function are obtained by imposing a scattering boundary 
condition. At the distance R where the nuclear potentials vanish, the upper component of 
the neutron radial wave function has the following asymptotic behavior: 



G v {kr) = A v \ji v (kr) - tan £,7^ (Ax)] , for r > R 



(24) 
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where ji u and n\ v are spherical Bessel and Neumann functions, respectively, and 5 U is the 
phase shift corresponding to the angular momentum (l„,j u ), k 2 = E 2 — M 2 . For the case 
of proton, the asymptotic behavior can be obtained by replacing the spherical Bessel and 
Neumann functions in Eq.(24) with the relativistic regular and irregular Coulomb wave 
functions 2j| , respectively. The energy of a resonant state is determined when the phase 
shift of the scattering state reaches tt/2. The wave function of scattering state is normalized 
to a delta function of energy 5(E — E'). This condition fixes the normalization constant A u , 
it can be obtained as: 

P llM + E , s 

A » = COs6 ^n-2J-> (25) 

In order to take into account the width effect of resonant states in the continuum we 
introduce a level density of the continuum into the gap equations. When we introduce 
single particle resonant states with widths into the pairing correlations, the gap equations 
can be expressed as: 



?( Ja+ 2) J(e a - A) 2 + A 2 + ^ ( j 



1W 9v{£v) de„=* (26) 



y/{e a - A) 2 + A 2 VJi* V(e u -\) 2 + A 2 G 



E u« + \ ) 1 - ^ ( j! A) 2 + A 2 +E ( ou + \ ) / 



y/(e v - A) 2 + A 2 



de v 
(27) 



The sums a and v run over the bound states and resonant states involved in the pairing 
calculation, respectively, and I„ is an energy interval associated with each partial wave 
il u ,3v)- The factor g v is defined as: 

9»(e v ) = , (28) 

7r de v 

which is the so-called continuum level density and 5 U is the phase shift of the scattering 
state with angular momentum v = (l u ,j u ). The factor g u represents the variation of the 
localization of scattering states in the energy region of a resonance, in other words, it reflects 
the widths effect of the resonant states. For a very narrow resonant state, the factor g v 
becomes a delta function. 

Taking account of those resonant continuum in the gap equations, the expressions of 
various densities in Eq.(16 - 19) have to be modified. They can be expressed as: 
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a 

2j„ + 1 



2 - y '" + 1 ^(|Ga(r)| 2 -|F a (r)| 2 ) 



47T 



(\G u (r)\ 2 -\F u (r)\ 2 )g u (e u )vl(e u )de„ , 



(29) 



^^ 2 (|G a (r)| 2 + |F a (r)| 2 ) 

2> + 1 
An 



(\G u (r)\ 2 + \F u (r)\ 2 )g u (e u )v 2 u (e„)de u , 



(30) 



p 3 (r) = x; 



2/ " ^ 2 (|G a (r)| 2 + |F a (r)| 2 ) 



47T 

a 

+ E^i^ r 3 / (|C,(r)| 2 + |F,(r)|V( £ > 2 ( £i/ )^, 



(31) 



2ja + l A-r 3 



47T 



^(|G„(r)| 2 +|F„(r)| 2 ) 



2j„ + 1 fl-T 3 

47T 



(|G,(r)| 2 + \F v {r)\ 2 )g v {e v )vl{e v )de v , (32) 



The Dirac equation Eq.(9) and meson field functions (12 - 14) as well as the photon 
field function (15) with corresponding densities (29 - 32) are solved self-consistently in an 
iterative way. Therefore, the total binding energy can be written as: 

e = J2 i 2 i" + !) v l E « + E ( 2 > + x ) / <?.(^K(^M^ 

-\ J {g*P s Vo + QuPuVo + g P PzP% + ep c A ) d 3 r 

E ( Ja + VaUa + E + y gv(£v)vv(s v )u v (e v )de v 



■-■41 -A' 1 / 3 
4 



(33) 



Where the last two terms are the pairing energy and the correction for the spurious center 
of mass motion, respectively. 
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IV. PROPERTIES OF NEUTRON RICH EVEN-EVEN NI ISOTOPES 

In this work we investigate the ground state properties of neutron rich even-even Ni 
isotopes in the RMF+BCS model and explicitly take account of the continuum. Partic- 
ularly we shall focus our attention on the effect of single particle resonant states on the 
pairing correlation in those neutron rich nuclei. Calculations are carried out in the RMF 
with the parameter set NL3. Three types of calculations: RMF+BCS, RMF+BCSR and 
RMF+BCSRW are performed. The effect of single particle states in the continuum on the 
ground state properties is investigated. 

In our investigation, the proton number of Ni isotopes is Z = 28, which is a closed shell. 
Therefore, the proton pairing gap is taken to be zero. For the neutron case, we use a state- 
independent pairing strength G = C/A, where the constant C = 20.5 MeV is adopted from 



Ref . 14] . It could also provide a best fit to reproduce experimental binding energies of known 
Ni-isotopes. In practical calculations, we restrict the pairing space to about one harmonic 
oscillator shell above and below the Fermi surface in the RMF+BCS model. The levels 
include l/ 5/2 , 2p 3/2 , 2p 1/2 ,lg 9/2 , 2d 5/2 , Ssy 2 , 2d 3/2 , lg 7 / 2 , lhn/ 2 and 2f 7/2 , some of them, 
such as 2d 5 / 2 , 2d 3 / 2 , lg7/ 2 , 2f 7 / 2 , and lh u / 2 are single particle resonant states, depending 
on the isotopes. Highly excited resonant states with large widths, such as lhg/ 2 and li\3/ 2 , 
are ignored in our calculations. Obviously, an enlarged pairing space in the normal BCS 
calculation would change the pairing contribution due to a constant pairing interaction. 
Therefore the results in cases of RMF+BCS and RMF+BCSR in Fig.l and Fig.2 would be 
different if one includes more states in the pairing active space. It is found that the pairing 
properties are not affected too much if one includes the resonant states lh 9 / 2 and l«i 3 / 2 with 
wide widths in the RMF+BCSRW calculation. 

The single particle energies and wave functions are first carried out by solving the Dirac 
equation Eq.(9) self-consistently. Using those single particle states we can solve the BCS 
gap equations in the three approaches. The Fermi energy and gap as well as the occupation 
probabilities of quasi-particle states are obtained simultaneously. The nuclear densities 
composed of quasi-particle states and potentials are recalculated. Then we solve Dirac 
equation again by an iterative way until the convergence is reached. The pairing correlation 
energies can be obtained in the Eq.(22) in the three approaches, which are shown in Fig.l 
for neutron rich even-even Ni isotopes of N = 40 to 70. Pronounced differences of the 
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pairing energies performed in various RMF+BCS approaches are observed for open shell 
nuclei, where the behavior of pairing energies is similar to those obtained by the RHB in 
Refs.[27]. The usual RMF+BCS approach produces the largest pairing energies, due to some 
non- resonant scattering states in the continuum are included. It can be seen in Fig.l that 
the pairing energies are reduced largely when the widths of single particle resonant states 
are taken into account in the RMF+BCSRW calculation. 

In order to understand the difference presented by different treatments on pairing cor- 
relations, we plot the results of pairing gaps and Fermi energies for neutron rich even-even 
Ni isotopes with masses from 68 to 98 in Fig.2(a) and Fig. 2(b). The curve in Fig. 2 (a) is 
the empirical pairing gap given by W.ljsfA. It is observed that pairing gaps are largely de- 
creased when the width effect of single particle resonant states is taken in account in pairing 
correlation calculations, while the Fermi energies remain unchanged for three a ppr oaches. 
This feature agrees with the case in the non-relativistic HF+BCS calculations [l0|. It is 
found that the pairing gaps in the RMF+BCSR calculation are close to those obtained in 
the box RMF+BCS approximation. As an example, for the neutron rich nucleus 90 Ni, the 
pairing gap and Fermi energy are A n =1.558 MeV and A n =-1.278 MeV in the RMF+BCSR 
approach, whereas, they are reduced to 1.223 MeV and -1.256 MeV, respectively when we 
perform the pairing correlation calculations by considering the width effect of resonances ex- 
plicitly. The reduced pairing gap and Fermi energy could change the occupation probability 
of neutrons at each single particle orbit near the Fermi surface. In Fig. 3 we show the occu- 
pation probabilities of neutron single particle levels in 90 Ni as a function of corresponding 
energies. The arrow is the position of the neutron Fermi surface. In 90 Ni the Fermi energy 
is very close to the continuum, therefore the continuum becomes important in the pairing 
correlations. Some scattering states in the continuum, such as states 3p3/2 and 3pi/2, which 
are not resonance states due to the low centrifugal barrier, are included in the RMF+BCS 
calculation. Therefore the RMF+BCS approach overestimates the pairing correlation and 
produces large pairing energies and pairing gaps. It is found that the width effect on the 
pairing is to reduce the pairing correlations, therefore occupation probabilities of low-lying 
states below the Fermi surface in the case of RMF+BCSRW are slightly larger than those 
without the width effect, and vice versa for states above Fermi surface. 



12 



In Fig. 4 we display the two-neutron separation energy S^: 

S 2n = B(Z,N)-B(Z,N-2) . (34) 

We calculate two-neutron separation energies for Ni isotopes till the neutron drip line in three 
approaches, RMF+BCS, RMF+BCSR and RMF+BCSRW. The two-neutron separation 
energies calculated in the RMF without pairings are also plotted in Fig. 4, which are denoted 
by crosses. The empirical data of S2n are obtained by Eq.(34), where the experimental 
binding energies of the Ni isotopes for iV <50 are taken from the Ref.j3|. The theoretical 
binding energies for Ni isotopes are listed in Table I. The position of the neutron drip line 
may be determined by the signature where the two-neutron separation energy changes its 
sign or the Fermi energy becomes positive. It is found that the position of the neutron drip 
line predicted by various RMF+BCS models is located between 98 Ni and 100 Ni, the same 
result is also obtained in RHB calculations[2j|. It is found that the two-neutron separation 
energies calculated in various RMF+BCS approaches are very close to each other, which 
can well reproduce the experimental data in the region of N = 42 to N = 50. In contrast, 
the results without pairing considerably deviate from the others. It indicates that the 
pairing is responsible to a reasonable description of the two-neutron separation energy. Some 
disagreements with experimental data are found for N = 38 and N = 40 isotopes, which 
were also observed in the RHB calculation with the parameter set NLsj^, Isflj] . Although 
pairing energies obtained in these three approaches with different treatment of pairings differ 
significantly, two-neutron separation energies are consistent with each other. This is because 
that differences, which appear in the pairing correlation energy are largely cancelled in the 

S2n- 

In table I, it is shown that the binding energies produced in the three RMF+BCS ap- 
proaches scatter when nuclei are far away from the j3 stable line. To clearly illustrate the 
width effect of the resonant states on the pairings we define a new quantity, which charac- 
terizes the pairing correlations energy: 

E BCS = E(RMF) - E(RMF + BCS) . (35) 

The Ebcs of Ni isotopes calculated in various approaches are plotted in Fig. 5. It is shown 
that the results produced in the three approaches are very close to each other for those 
isotopes not far away from the j3 stable-line. It implies that the width of the resonance in 
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the pairing correlation could be neglected at iV < 34. The RMF+BCS and RMF+BCSR give 
more or less similar values of E BC s, even near the drip line. The width effect gets more and 
more pronounced as the neutron number increases, especially near the drip line. Therefore 
a proper treatment of the resonant continuum including its width might be necessary for 
the nuclei near the drip line. 

Neutron rms radii of neutron rich even-even Ni isotopes are further calculated in the 
three approaches, which are plotted in Fig. 6. The rms radii calculated in the RMF are also 
shown in Fig. 6 and displayed by cross. It is seen that neutron rms radii of Ni isotopes at 
A < 80 calculated in all approaches are very close to each other, which is similar to what was 
observed in the binding energy. The neutron rms radii calculated in the RMF+BCS with a 
box become larger at A > 80 in comparison with those produced in other approaches. This 
is due to the fact that the contribution of the pairing correlation from the scattering states 
3p 3 / 2 and 3pi/ 2 is included in the RMF+BCS calculations. They are not resonant states 
and their wave functions are not well localized inside the nucleus. The neutron rms radius 
for 84 Ni given in the RMF is slightly smaller than that obtained in the other approaches. 
Actually this is a natural result when the pairing correlation is switched on. The neutron 
rms radius in the RMF+BCS is increased because the pairing interaction scatters some 
neutrons from 2d 5 / 2 to the loosely bound state Ss 1 / 2 which is a state less localized. In 
contrast, in the case of 90 Ni the states 3si/ 2 and 2d 3 / 2 , which are completely occupied in the 
RMF approximation are scattered to states 1*77/2 and lhn/2 due to the pairing interaction. 
Although the states lc/7/2 and l/in/2 are close to or buried in the continuum, their wave 
functions are more localized inside the nucleus than those for the 3si/ 2 and 2d 3 / 2 states due 
to a high centrifugal barrier. Therefore the rms radius of 90 Ni calculated in the RMF is 
larger than that produced in calculations including the pairing correlation. 

We plot the neutron density of 84 Ni in Fig. 7, where the dash-dotted curve is obtained 
in the RMF, the short-dotted, dashed and solid curves are produced in the RMF+BCS, 
RMF+BCSR and RMF+BCSRW approaches, respectively. The tail of the density distribu- 
tion gets larger when the pairing correlation is considered. The width effect on the density 
distribution is very small, which is consistent to the calculated neutron rms radii. In order 
to show the effect of the box size on the density in the box RMF+BCS approximation we 
also plot the neutron density with various values of the box size Rb ox = 15, 20, 25 fm in 
Fig. 7. It is observed that the tail of the density strongly depends on the box size, where 
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unphysical particle gas may appear in the exotic nucleus. In the RMF+BCS calculation 
non-resonant discretized states in the continuum are included in pairing correlation, such 
as 3p3/2 and 3pi/2 states. Their wave functions have a long tail strongly depending on the 
box size. It is found that this problem is well overcome when one performs the pairing 
correlation calculation with only a few narrow resonant states in stead of those discretized 
states in the continuum. 



V. SUMMARY 



In this paper, we have investigated the pairing correlation for neutron-rich Ni isotopes 
in the relativistic mean field theory and Bardeen-Cooper-Schrieffer approximation with a 
constant pairing strength. A proper treatment of the single particle resonant state in the 
continuum on pairing correlations has to include not only its energy, but also its width. 
The inclusion of the width of the resonant state in the pairings can be performed by in- 
troducing a level density in the continuum into the pairing gap equation. The resonant 
continuum is solved by imposing a proper scattering boundary condition. The investigation 
is performed in three approaches: RMF+BCS, RMF+BCSR and RMF+BCSRW with the 
effective Lagrangian parameter set NL3. A special attention is paid on the width effect 
of resonant states in the continuum on the pairing correlation for nuclei close to the drip 
line. We have studied the width effect of the resonant continuum on pairings: pairing gap, 
Fermi energy and occupation probability, as well as nuclear ground state properties, such 
as binding energy, two-neutron separation energy, neutron rms radius and neutron density. 
The results show that a proper treatment of the resonant continuum in pairing correlations, 
which explicitly includes the width of the single particle resonant state is important for 
nuclei close to neutron drip line. They could affect the pairing gaps, Fermi energies, pairing 
energies, and binding energies considerably. Various RMF+BCS approaches could give a 
similar description on the ground state properties for nuclei not far away from the (3 stable 
line. It is observed that unphysical particle gas, which may appear in the traditional mean 
field plus BCS calculation for nuclei in the vicinity of drip line can be well overcome when 
one performs pairing correlation calculations with only resonant states instead of discretized 
states in the continuum. It might be concluded that the simple BCS method works well in 
describing the pairing correlations for neutron rich nuclei provided the continuum is properly 
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treated. 
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TABLE I: The binding energies of Ni isotopes near the neutron drip line calculated in the 
RMF+BCS, RMF+BCSR, and RMF+BCSRW approaches with parameter set NL3. The available 
experimental binding energies of the Ni isotopes for N <50 from Ref. 
values in the table are in unit of MeV. 



28 1 are also shown. All energy 





RMF+BCS 


RMF+BCSR 


RMF+BCSRW 


Exp. 


58 M 


502.347 


502.275 


502.104 


506.453 


60 M 


521.655 


521.448 


521.363 


526.841 


62 Ni 


540.458 


540.374 


540.251 


545.258 


64 Ni 


558.361 


558.273 


558.055 


561.754 


66 Ni 


575.642 


575.579 


575.009 


576.830 


68 Ni 


590.892 


590.960 


590.930 


590.430 


70 Ni 


603.233 


603.195 


602.775 


602.600 


72 Ni 


614.246 


614.116 


613.567 


613.900 


74 Ni 


624.428 


624.307 


623.752 


623.900 


76 Ni 


633.840 


633.816 


633.414 


633.100 


78 Ni 


642.352 


642.555 


642.555 


641.400 


80 M 


647.837 


647.659 


646.964 




82 M 


652.279 


651.943 


650.888 




84 M 


656.092 


655.678 


654.472 




86 M 


659.422 


658.975 


657.682 




88 M 


662.284 


661.863 


660.582 




90 Ni 


664.756 


664.319 


663.047 




92 Ni 


667.756 


666.456 


665.254 




94 Ni 


668.562 


668.319 


667.144 




96 Ni 


670.133 


670.183 


669.332 




98 M 


671.382 


671.605 


671.605 
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FIG. 1: The pairing energies for neutron rich even-even Ni isotopes with N = 40 — 70. The solid 
squares, circles and left-triangles denote the results calculated in the RMF+BCS, RMF+BCSR 
and RMF+BCSRW approaches, respectively. All calculations are performed with the parameter 
set NL3. 
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FIG. 2: The neutron pairing gap (a) and Fermi energy (b) as a function of the atomic number A. 
The notations are same as in Fig.l. 
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FIG. 3: The occupation probabilities for Ni as a function of the single particle energy around the 
threshold. The results are performed in RMF+BCS, RMF+BCSR and RMF+BCSRW approaches 
with the effective nonlinear interaction NL3. The notations are same as in Fig.l. 
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FIG. 4: Two-neutron separation energies for neutron rich even-even Ni isotopes calculated in the 
RMF, RMF+BCS, RMF+BCSR, and RMF+BCSRW approaches with parameter set NL3. The 
cross signs are the results calculated in the RMF. The experimental data denoted by open squares 
are taken from Ref. 28 1. Other notations are same as in Fig.l. 
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FIG. 5: The pairing correlation energies Ebcs calculated in Eq.(35) in the RMF+BCS, 
RMF+BCSR and RMF+BCSRW approaches. The notations are same as in Fig.l. 
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FIG. 6: Neutron rms radii for neutron rich even-even Ni isotopes calculated in the RMF, 
RMF+BCS, RMF+BCSR, and RMF+BCSRW approaches with parameter set NL3. The notations 
are same as in Fig.4 
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FIG. 7: Neutron density distribution for 84 Ni calculated in the RMF, RMF+BCS, RMF+BCSR, 
and RMF+BCSRW approaches, where the RMF+BCS results calculated with various values of 
the box size Rb ox = 15, 20, 25 fm are also plotted. 



